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Abstract: Utilizing only the vertical muon intensity of the Gaisser parameterization of the muon flux 
at the surface and propagating this energy spectrum underground according to statistical ionization and 
radiative energy losses, it is possible to calculate the underground muon intensity Crouch curve. In addi- 
tion, the primary spectral index of the Gaisser parameterization can be adjusted from E~ to E~ 
simply by minimizing the deviation from the Crouch curve. For chemical compositions other than stan- 
dard rock, the propagation of the spectrum underground can be repeated with a different muon energy 
loss in the material. The resulting underground muon intensity curve represents a consistent conversion 
of the Crouch curve to the local rock, fully accounting for the energy dependence of the muon dE/dx. 



Gaisser Parameterization of Muon Flux 
at Surface 

The Gaisser parameterization of the cosmic in- 
duced muon flux at surface is an approximate ex- 
trapolation formula valid when muon decay is neg- 
ligible {Ef^ > 100/ cos 6* GeV) and the under- 
ground detector zenith angle 6 can be assumed 
identical to the production angle in the upper at- 
mosphere {6 < 70°). Not including small contri- 
butions from charm and heavier flavors, which are 
negligible for energies below lOTeV, the Gaisser 
parameterization is [1][2], 
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where the two terms give the contribution of pions 
and charged kaons. The energy spectrum steep- 
ens gradually to reflect the primary spectrum in 
the 10 — 100 GeV range, and steepens further at 
higher energies because pions with Et, > critical 
energy e^r ~ 115 GeV tend to interact in the atmo- 
sphere before they decay. Asymptotically (E^ ^ 
1 TeV), the energy spectrum of atmospheric muons 
is one power steeper (E^^ '^) than the primary 



spectrum (i?^^-^). For the following calculations 
we only consider the vertical muon intensity for 
which cos 9 = 1. 



Muon Range Underground 

The statistical energy loss of muons, traversing an 
amount X of matter, with energies far above the 
Bethe-Bloch minimum is given as 
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where a is the coUisional term (i.e. ionization 
mostly due to delta-ray production) and the second 
term accounts for the three radiative muon energy 
loss processes: 1 . Bremsstrahlung, 2. pair produc- 
tion and 3. photonuclear interactions. In Table 1 
these energy loss parameters are listed for standard 
rock. The critical energy where ionization losses 
equal radiative losses in standard rock is approxi- 
mately 0.6 TeV. As the energy dependencies of the 
a and b parameters are relative mild, they are often 
assumed to be constant and the then resulting sim- 
plified differential equation can be easily solved by 
using an exponential function [2]. However, for 
our purpose of precisely determining the average 
muon range underground for each value of surface 
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energy, the a and S6 values are parameterized with 
the following functions 



Ao + iA,-\og,oE^[GeV]) (3) 
Bo + (Bi Aog^o E^[GeV]) + 
{B2-i\og,oE^[GeV]f}. (4) 



Table 2 lists the fitted coefficients for standard 
rock. In addition, a Geant4 [4] based simulation 
is performed, tracking muons through a block of 
standard rock, defined as a mixture of CaCOs and 
MgCOz (with a density of 2.65 gjcrn? and mass 
fractions of 52% O, 27% Ca, 12% C and 9%Mg). 
The comparison for the simulated parameter as 
a function of energy and the parameterization in 
Eq. 4 validates that the fractional difference is less 
than 3 % over the entire energy range. 
We used Eq. 2 to numerically compute the propa- 
gation of the muon by stepping (AX = 1 g/cm?) 
through standard rock. Thus, for each initial value 
of muon energy, we determined, at what value of 
depth in meter-water-equivalent the muon ranges 
out. Fig. 1 shows the computed muon range in 
mwe of standard rock as a function of muon sur- 
face energy from lOGeV to lOTeV. Finally, the 
muon range X in standard rock can be simply pa- 
rameterized as 



X[mwe] = po ■ logJ(pi • E^[GeV]) +P2}, (5) 

with the fitted parameters listed in Table 3. 
The same computation can be done for overbur- 
dens made up of other rock compositions like e.g. 
for more dense rock (p « 2.85 g/cm^) at the 
Soudan iron mine in northern Minnesota, US (c.f. 
again Fig. 1 and Table 3). In order to easily trans- 
form the muon energy loss parameters and bg in 
standard rock (average nuclear properties = 11, 
As = 22 [2]) to a' and b' in other rock (e.g. Soudan 
= 12, ]?' = 24 [5]), the a and total b parameter 
at a given muon energy can be scaled as 
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The muon energy loss from Bremsstrahlung and 
pair production can be scaled with Z'^ /A to the 
first order [5], whereas the simulated contribution 
from photonuclear interactions to dE / dX has only 
a very weak dependence on the nuclear properties 
in rock and is herein assumed constant per mwe. 
Furthermore, the MC simulation shows that in the 
region of interest from 250 GeV to 10 TeV the pho- 
tonuclear interactions account for a constant frac- 
tion (10 %) of all radiative muon energy losses in 
standard rock. The scaled a and S6 values can be 
parameterized again with the functions in Eq. 3, 4 
and the average muon range underground for each 
value of surface energy can be computed again in 
the same procedure as described above {c.f. again 
Table 2 and 3). 
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Table 1 : Average muon energy loss parameters calcu- 
lated for standard rock [3] [2] 
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Table 2; Fitted coefficients for the parameterizations of 
the ionization (Eq. 3) and the total radiative (Eq. 4) 
average muon energy loss in standard and Soudan rock. 
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Table 3: Fitted coefficients for the parameterization 
(Eq. 5) of the average muon range in standard and 
Soudan rock. 



Calculating the Muon Intensity Under- 
ground 

Together with Eq. 1, which is the Gaisser pa- 
rameterization of the differential intensity of ver- 
tical muons at the surface (cos^ = 1), and Eq. 
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Figure 1: Calculated average muon range in standard 
and Soudan rock as a function of initial muon energy at 
the surface. 

5, which relates the average muon range under- 
ground to a given value of surface energy, it is 
now very simple to compute the intensity under- 
ground as a function of slant depth in standard 
rock. For this, the differential intensity of verti- 
cal muons at the surface (Eq. 1 with cos 9 — 1) is 
stepwise added (AE^ = lOMeV) from 100 TeV 
down to 1 GeV. At each step the average muon 
range X in mwe is calculated with Eq. 5 and as- 
signed to the interim value of the intensity sum. 
Thus, the integral vertical muon intensity under- 
ground as a function of slant depth in standard rock 
can be efficiently computed from 10000 mwe up 
to the surface. Fig. 2 depicts the result in compar- 
ison with the Crouch curve, which is the param- 
eterized 'world average' of deeply underground 
measured (> 1000 mwe) vertical muon intensi- 
ties; compiled by Crouch and referring to standard 
rock [6]. The result of our calculation for stan- 
dard rock agrees already well with the underground 
measured Crouch curve. Both the intensity values 
and the spectral shape of the exponentially falling 
off Crouch function are well reproduced. At the 
small slant depth value of 1000 mwe, where under- 
ground measurements start to be considered in the 
Crouch fit, both curves match very well, whereas 
for larger values of slant depth the new calculation 
tends to fall off slightly faster. 



Figure 2: Calculated integral vertical muon intensity 
underground as a function of slant depth in standard rock 
in comparison with the Crouch function [6]. 

Refining the Gaisser Parameterization 
by Comparison with Crouch Curve 

We tried to further improve the agreement between 
the calculation and the Crouch curve by refining 
the primary spectral index (7) and the absolute nor- 
malization (c) of the Gaisser parameterization of 
the muon flux at the surface (substitute i?^^-^ with 
c ■ E'' in Eq. 1). In a Newtonian iteration proce- 
dure, the two variables 7 and c were varied before 
each new calculation, such that the new values at 
^ 1500 mwe and ^ 9500 mwe are in optimum 
agreement with Crouch and the amplitude of the 
fractional differences for slant depth values in be- 
tween is minimal (c.f. Fig. 3). The best fit values 
obtained by this procedure are 7/^4 — —2.643 for 
the primary spectral index and Cnorm — 80.5% 
for the normalization factor (c.f. Fig. 4, 2). We 
found the maximum fractional difference between 
the optimized calculation and the Crouch curve to 
be 6 %, which is of the order of the uncertainty on 
the measured Crouch curve. 

Discussion 

The fact that the parameterized muon flux at the 
surface in combination with the propagated muon 
energy loss in the rock reproduces the 'world av- 
erage' Crouch curve implies that the underlying 
physics, as described in this paper, is well un- 
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Figure 3: Fractional difference between tiie integral 
vertical muon intensity underground according to the 
Crouch curve and the optimized calculation for slant 
depth values in standard rock. 

derstood. Our optimized Gaisser parameteriza- 
tion, with an absolute value of the primary spec- 
tral index of slightly less than 2.7 and a ^ 20% 
lower muon intensity at surface, is in reasonable 
agreement with experimental data [2]. Crouch [6] 
converted the underground muon intensity data to 
depth values in standard rock utilizing an approxi- 
mate mathematical solution of the differential Eq. 
2 assuming the a and b parameters to be constant 
[2], which becomes less valid at increasing depths. 
It is possible now to perform the herein described 
calculation in other rock than standard rock, yield- 
ing new underground muon intensity curves for the 
local rock composition. 

Summary 

Taking the vertical muon intensity of the refined 
Gaisser parameterization of the muon flux at the 
surface (80.5% • i;-2-643 instead of E-^-'^) and 
propagating this energy spectrum underground ac- 
cording to statistical ionization and radiative en- 
ergy losses yields a good fit to the underground 
muon intensity Crouch curve. The obtained agree- 
ment is for most part better than the uncertainty of 
~ 5 % associated with the Crouch function. For 
chemical compositions other than standard rock, a 
consistent computation of the underground muon 
intensity curve can be repeated, fully accounting 



Figure 4; Integral vertical muon intensity underground 
according to the Crouch curve (blue) and the optimized 
calculations for slant depth values in standard (green) 
and Soudan rock (black), as well as an approximate con- 
version (red) of the Crouch curve to Soudan rock (con- 
stant muon energy loss parameters a and 6 at 1 TeV). 

for the energy dependence of the muon dE/dx in 
the local rock. This can yield a better determina- 
tion of a map of the overburden of an underground 
muon detector, by normalizing the measured in- 
tensity to the computed distribution as a function 
of slant depth in the local rock. Furthermore, the 
methods described herein also allow for a more 
precise extrapolation of the underground muon en- 
ergy back to the surface. 
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